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1.  Proposed  summary 

The  project  was  designed  to  develop  novel  spectroscopic  tools  for  the  determinations  of 
molecular  structures  and  dynamics  on  the  surfaces  of  metal  nanomaterials  -  the  critical 
component  of  heterogeneous  catalysts,  and  to  use  the  new  tools  to  investigate  whether  a 
fundamental  physical  principle,  the  Bom-Oppenheimer  Approximation  (BOA),  is  valid  on  the 
surfaces  of  these  materials  or  not. 

2.  Achievements 

During  the  past  three  years  supported  by  the  award  (. FA9550-11-1-0070 ),  my  group 
pioneered  two  novel  molecular  structural  tools:  (1)  the  vibrational  cross  angle  method  for  the 
determination  of  3D  molecular  structures1'5;  and  (2)  the  vibrational  energy  transfer  method  for 
the  determination  of  intermolecular  distances  with  an  angstrom  precision6'11.  Using  the  methods, 
we  showed  that  the  3D  molecular  conformations  on  metal  nanomaterials  with  Au  nanoparticles 
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as  examples  can  be  directly  determined  ’  ,  and  the  validity  of  Bom-Oppenheimer  Approximation 
on  the  surfaces  of  metal  nanomaterials  is  dependent  on  detailed  molecular  properties5’12'14,  and 
the  energy  dissipation  pathways  in  the  systems  largely  depend  on  whether  BOA  is  valid12'14.  In 
the  following,  the  accomplishments  will  be  briefly  described. 

2.1  The  vibrational  cross  angle  method  to  determine  3D  molecular  structures 

One  of  the  major  problems  in  experimentally  studying  heterogeneous  catalysis  is  the  lack 
of  tools  for  the  determinations  of  molecular  structures  on  the  surfaces  of  metal  nanomaterials. 
Practical  catalysts,  e.g.  oxide-supported  metal  clusters,  are  polycrystalline  solids.  XRD  cannot  be 
used  to  resolve  the  molecular  structures  in  such  systems  because  they  are  not  single  crystals. 
NMR  also  has  difficulties  in  determining  molecular  structures  in  them  not  only  because  the  solid 
NMR  has  broad  lineshapes,  but  also  because  the  metal  clusters  typically  have  open  shell 
electrons  which  can  severely  disrupt  the  magnetic  fields  on  their  surfaces.  ID  vibrational 
methods,  e.g.  FTIR,  Raman,  or  SFG,  are  usually  used  in  monitoring  the  molecular  changes  in 
heterogeneous  catalytic  systems  by  obtaining  the  vibrational  frequencies  and  peak  intensities. 
However,  the  structural  determination  capability  of  the  ID  vibrational  techniques  suffers  from 
the  fact  that  vibrational  frequencies  are  severely  affected  by  Fermi  resonances  because  of  which 
the  vibrational  frequency  shifts  can  be  more  significant  than  those  caused  by  molecular  structural 
changes.  Motivated  by  such  a  practical  difficulty,  my  groups  developed  a  vibrational  cross  angle 
method  which  is  general  for  the  determination  of  3D  molecular  structures  in  condensed  phases. 

The  principle  of  our  method  is  by  determining  the  cross  angles  between  vibrational 
transition  dipoles  covering  the  entire  molecular  space  to  determine  3D  molecular  conformations 
(Fermi  resonances  don’t  affect  the  results  from  such  a  method),1'5  with  the  aid  of  our  unique 
powerful  ultrafast  multiple  dimensional  IR/THz  spectroscopy  setup3  which  solves  a  longstanding 
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problem:  the  lack  of  frequency  tunable  sources  over  the  entire  Mid-IR  range  with  sufficient 
power  has  prevented  essentially  all  other  current  2D  IR  techniques  from  resolving  3D  molecular 
structures  because  vibrational  couplings  with  wide  spatial  or  frequency  separation  cannot  be 
observed  or  modes  with  weak  transition  dipole  moments  are  too  weak  to  be  detected.  The 
method  has  been  benchmarked  with  well-defined  crystalline  samples,1’4  and  successfully  applied 
to  map  the  3D  molecular  conformations  in  liquids,1  solids1,4  and  on  the  surfaces  of 
nanomaterials5. 

Our  recent  work1  will  be  used  to  demonstrate  the  method. 


Figure  1.  Molecular  conformations  ofMNA  in  (A)  the  white  crystal  MNA-W,  and  (B)  the  yellow 
crystal.  (C)  Illustration  of  how  the  vibrational  cross  angle  between  two  modes  is  experimentally 
determined. 

Fig.l(A)&(B)  display  the  molecular  conformations  of  MNA  (4'-Methyl-2'- 
nitroacetanilide)  in  two  of  its  crystals  (white  and  yellow  as  shown)  determined  by  XRD.  The 
conformations  are  defined  by  the  two  dihedral  angles  ACC  /  NO  (between  the  benzene  plane 
and  the  nitro  plane)  and  ACC  INC  (between  the  benzene  plane  and  the  amide  plane).  In  the 
white  crystal,  there  are  one  pair  of  chiral  images,  and  in  the  yellow  crystal  there  are  two  pairs. 
The  molecule  has  many  typical  chemical  bonds,  e.g.  C  =  C-H ,  N  =  0 ,  C  =  C,  C  =  0 ,  N-H 
and  C-0 ,  which  cover  the  three  planes  of  the  molecule.  The  relative  orientations  of  the 
vibrations  of  these  bonds  can  be  utilized  to  obtain  the  molecular  conformations  in  terms  of  the 
two  dihedral  angle  values. 

To  determine  the  cross  angle  between  two  vibrations,  as  illustrated  in  fig.  1(C),  a  linearly 
polarized  IR  pulse  excites  a  vibrational  mode  (the  NO2  stretch  at  1362  cm'1).  After  a  very  short 
period  of  time  (0.1~0.2ps)  before  the  molecular  rotation  or  conformational  changes  have 
occurred  for  a  substantial  extent,  another  linearly  polarized  pulse  of  different  frequencies  detects 
a  signal  generated  from  the  response  of  another  vibrational  mode  (the  C=0  stretch)  to  the 
excitation  of  the  NO2  stretch.  In  general,  the  excitation  of  one  vibrational  mode  can  lead  to  the 
vibrational  frequency  shift  of  another  mode  because  of  the  anharmonic  coupling.  The  coupling 
produces  a  cross  peak  pair  in  the  experimental  results  as  displayed  in  fig.2B.  In  fig.2B,  the 
excitation  frequency  cox  =1362  cm"1  is  the  NO2  stretch  0-1  transition  frequency,  indicating  that 
the  cross  peak  pair  is  from  the  NO2  excitation.  The  detection  frequency  (red  peak) 
a>3  =  1672  cm”1  is  the  C=0  stretch  0-1  transition  frequency,  indicating  that  the  cross  peak  pair  is 

from  the  0-1  transition  frequency  shift  of  the  C=0  stretch  from  1672  cm'1  (red  peak)  to  1555  cm' 
1  (blue  peak)  caused  by  the  NO2  excitation.  The  amplitudes  of  the  cross  peaks  (fig.2D)  are 
dependent  on  the  polarizations  of  the  exciting  and  detecting  beams,  and  the  cross  angle  6 
(0°  <6  <  90°)  between  the  transition  dipole  moment  directions  of  the  two  coupled  modes.  For  a 
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sample  isotropically  distributed  within  the  laser  focus  spot,  the  vibration  cross  angle  can  be 
straightforwardly  determined  based  on  the  relation: 

/,  2-cos2#  . 

—  = - eq.l 

I{1  1  +  2  cos2# 

where  In,I±  are  off-diagonal  peak  intensities  with  parallel  and  perpendicular  excitation/detection 
polarizations,  respectively.  From  both  parallel  and  perpendicular  measurements  (fig.2D)  and  eq.l, 
the  vibrational  cross  angle  between  the  NO2  stretch  and  the  C=0  stretch  is  determined  to  be 
70°  ±2°. 


Figure  2.  (A)  Multiple-mode  2D  IR  spectrum  of  a  polycrystalline  sample  of  the  MNA  white 
crystal.  (B)  Enlarged  2D-IR  spectrum  for  the  cross  peak  pair  N02/C=0,  and  (C)  enlarged  2D-IR 
spectrum  for  the  cross  peak  pair  NH/C=0.  (D)  A  slice  cut  along  cox  =1362  cm1  of  fig.2B  with 
the  polarization  of  the  excitation  both  parallel  (\\)  and  perpendicular  (Z)  to  the  polarization  of 
the  detection  beam  and  (E)  A  slice  cut  along  cox  =3260  cm1  of  fig.2C  with  both  parallel  and 
perpendicular  configurations.  The  average  difference  between  the  experimental  and  calculated 
vibrational  cross  angles  of  MNA  conformations  with  different  ZCC  /  NO  and  ZCC  /  NC 
dihedral  angles  for  (F)  MNA-W,  (G)  MNA-Y1,  and  (H)  MNA-Y2.  The  z-axis  is  the  amplitude  of 
difference.  The  minimum  difference  values  are  labeled  with  white  boxes.  The  dihedral  angles  of 
red  dots  are  determined  by  XRD. 

Through  a  similar  procedure  by  scanning  frequencies  from  1200  cm'1  to  3400  cm'1 
(fig.2A),  the  cross  angles  among  other  modes  were  also  obtained.  For  each  sample,  we  collected 
vibrational  cross  angles  of  12~15  pairs  of  coupled  modes.  These  experimentally  obtained 
vibrational  cross  angles  are  different  from  the  chemical  bond  cross  angles  needed  for 
constructing  3D  molecular  conformations.  Theoretical  calculations  are  then  used  to  convert  the 
vibrational  angles  into  bond  angles.  In  doing  so,  we  preset  a  series  of  molecular  conformations 
and  calculated  their  vibrational  cross  angles  and  compared  the  calculated  angles  to  the 
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experimentally  measured  angles  until  a  minimum  was  found.  The  conformations  with  calculated 
vibrational  angles  closest  to  the  measured  values  were  then  considered  as  the  most  probable 
conformations  experimentally  determined.  The  difference  between  the  calculated  and  measured 
vibrational  cross  angles  versus  the  two  dihedral  angles  is  plotted  in  fig.2(F)~(H)  for  the  three 
pairs  of  chiral  images  in  the  two  crystals  (The  two  pairs  of  chiral  images  in  the  yellow  crystal  can 
be  individually  determined  because  they  have  distinct  vibrational  frequencies  for  some 
vibrational  modes).  The  dihedral  angles  (marked  with  white  boxes  in  fig.2(F)~(H))  determined  in 
this  way  are  very  similar  to  those  determined  by  XRD  (marked  with  red  dots).  For  the  same 
molecule,  we  also  measured  its  conformations  in  liquids  with  a  similar  approach.  The  determined 
liquid  conformations  are  very  similar  to  those  from  theoretical  calculations.  Using  the  method, 
the  conformation  of  a  ligand  molecule  on  the  surface  of  4nm  Au  particle  was  also  determined 

<fig.3).5 


Figure  3.  Illustration  of  molecular  structure  on  the  surface  of  4nm  Au  particle  is  determined 
with  the  vibrational  cross  angle  method. 

2.2  Vibrational  energy  transfer  method  to  determine  angstrom  molecular  distances 

The  nature  of  the  Bom-Oppenheimer  Approximation  (BOA)  is  that  electrons  are  much 
lighter  than  nuclei  so  that  electronic  motions  are  assumed  to  be  much  faster  than  nuclear  motions. 
In  our  experiments,  the  criterion  to  judge  the  validity  of  BOA  is  whether  electronic  and  nuclear 
motions  are  coupled.  If  they  are  strongly  coupled,  the  vibrational  motions  prepared  by  laser 
excitation  can  be  converted  into  electronic  motions  very  fast,  resulting  in  experimentally 
observed  fast  vibrational  relaxations.  However,  vibrational  relaxations  in  condensed  phases  can 
also  be  induced  by  vibration/vibration  couplings.  To  unambiguously  distinguish  these  two 
origins,  we  need  to  understand  how  vibrational  energy  transfers  because  of  vibration/vibration 
couplings.  Historically,  vibrational  energy  transfers  or  relaxations  induced  by  vibration/vibration 
couplings  were  described  with  the  correlation  formalism  or  the  phonon  compensation 
mechanism  or  the  dark  mode  formalism.  However,  none  of  the  previous  theoretical  approaches  is 
experimentally  useful,  as  they  contain  too  many  parameters  that  cannot  be  experimentally 
determined.  Therefore,  the  first  step  of  our  experiments  is  to  design  experiments  to  search  for  the 
quantitative  description  of  vibrational  energy  transfers  in  condensed  phases  with  all  parameters 
that  can  be  experimentally  determined.  From  a  series  of  systematic  studies,  we  have  developed  a 
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general  theory  of  which  the  Forster  theory  is  only  an  extreme  case  to  quantitatively  describe 
intermolecular  vibrational  energy  transfers6,7.  The  research  led  to  the  birth  of  the  vibrational 
energy  transfer  angstrom  molecular  rule  method  ’  .  The  method  was  benchmarked  with  a 
crystalline  sample  with  known  well-defined  donor/acceptor  distances,  as  shown  in  fig.46’7,  and 
successfully  applied  to  investigate  ion  clustering  in  strong  electrolyte  aqueous  solutions.15,16  The 
results  demonstrated  that  cations  directly  bind  to  anions,9  that  mixed  ions  can  segregate  because 
of  different  water  affinities,10’15  and  that  organic  molecules,  that  are  models  of  protein  building 
blocks,  can  directly  bind  to  strong  denaturant  anions.11  In  contrast  to  electronic  “FRET”  which 
uses  large  chromophores  that  can  distort  structure,  the  vibrational  energy  transfer  method  can  use 
labels  that  have  no  effect  upon  structure.  Its  chromophores  are  the  chemical  bonds  themselves. 
From  our  research  ’  ,  we  also  found  that  the  previous  phonon  compensation  mechanism  to 
describe  intermolecular  vibrational  energy  transfers  in  liquids  was  wrong.  The  dominant 
vibrational  energy  transfer  mechanism  in  liquid  should  be  our  dephasing  mechanism.  This  is 
because  the  spectral  diffusion  is  much  faster  than  the  energy  transfers  so  that  the  modulations  by 
phonons  (or  instantaneous  normal  modes)  on  both  energy  donor  and  acceptor  cancel  out.  In  other 
words,  phonons  have  no  impact  on  the  energy  transfer  process.  The  energy  transfer  driving  force 
is  the  frequency  fluctuation  because  of  the  dephasing  events  caused  by  molecular  collisions.  In 
solids,  phonon  modes  are  long-lived,  the  dephasing  mechanism  and  the  phonon  compensation 
mechanism  are  competing  with  each  other.  Even  for  the  transfers  in  solids,  we  found  that  the 
majority  of  previous  researchers  made  a  mistake  in  taking  the  first  order  coupling  matrix  to 
describe  the  phonon  compensation  mechanism  for  the  one  phonon  process.  Instead,  the  term 
playing  a  significant  role  should  be  the  second  order  coupling  matrix.  With  our  theory,  the 
experimentally  observed  surprising  opposite  temperature  dependences  of  resonant  and 
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Figure  4.  Illustration  of  determining  the  shortest  distance  between  two  anions  in  KSCN  crystal 
with  the  vibrational  energy  transfer  method. 

2.3  Electron/vibrational  coupling  on  surfaces  of  metal  nanomaterials 

Equipped  with  the  new  spectroscopic  setup  and  the  knowledge  about  vibrational  energy 
transfers,  we  were  ready  to  investigate  the  BOA  (or  electron/vibrational  coupling)  problem  on 
metal  nanomaterial  surfaces12'14.  First  of  all,  we  found  that  the  electron/vibrational  couplings  are 
particle  size  dependent  (fig.5):  when  the  particle  is  equal  to  or  larger  than  2nm,  a  strong 
electron/vibrational  coupling  is  observed.  When  the  particle  is  smaller  than  2nm,  no  clear 
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electron/vibrational  coupling  is  observed.  The  size  dependence  is  the  same  for  Pt,  Pd,  and  Au 
particles.  However,  the  electron/coupling  strengths  are  very  similar  in  Pt  and  Pd  both  of  which 
are  much  larger  than  that  of  Au  (fig.6).  The  observation  can  be  fully  explained  by  the  facts:  (1) 
particles  larger  than  2nm  are  already  metallic  so  that  their  surface  electronic  levels  are 
continuous  which  can  be  on  resonant  with  the  surface  molecular  vibration  (CO  stretch  in  our 
studies).  Particles  smaller  than  2nm  are  semiconductive  and  their  electronic  energy  levels  are 
much  higher  than  that  of  the  surface  molecular  vibrations  so  that  resonant  energy  transfers 
between  surface  electrons  and  vibrations  cannot  occur.  (2)  Pt  and  Pd  have  similar  interaction 
strengths  with  surface  molecules  (CO)  by  donating  d  electrons  into  the  antiorbitals  of  the  CO 
molecules,  which  is  much  stronger  than  that  between  Au  and  CO  molecules  in  which  the  d 
electron  donation  doesn’t  exist.  The  electron/vibrational  coupling  is  also  surface  site  dependent. 
On  2nm  pt  particle  surface  (fig.7),  the  order  of  coupling  strengths  is  hollow>bridge>step>  terrace. 
A  similar  site  dependence  is  also  observed  on  Pd  surfaces.  Surface  distance  dependent 
measurements  show  that  the  energy  relaxation  rate  is  approximately  proportional  to  the  2.5th 
power  of  the  distance  between  the  vibrational  mode  and  particle  surface,  indicating  that  the 
electron/vibrational  coupling  induced  relaxation  is  mainly  through  the  through-bond  mechanism 
(the  “mechanical”  mechanism  rather  than  the  “dipole/dipole”  mechanism).  Temperature 
dependent  measurements  from  80-3 00K  show  that  the  relaxation  induced  by  electron/vibrational 
coupling  is  temperature  independent,  strongly  suggesting  that  this  is  a  resonant  transfer  process. 
Because  of  the  size  dependent  electron/vibrational  coupling  strengths,  the  heat  generation 
pathways  on  the  particles  are  also  size  dependent.  On  Pt  particles  larger  than  2nm,  heat  is 
generated  within  a  few  ps,  while  it  is  longer  than  50ps  on  the  lnm  particle  (fig.8).  Based  on  the 
results,  the  energy  dissipation  pathways  are  summarized  in  fig.9. 


Figure  5.  The  relaxations  of  CO  1st  excited  state  excitation  (red:  0-1  transition;  blue:  1-2 
transition)  on  the  surfaces  of  Pt  particles  (step  surface  site)  with  different  sizes:  (A)  lnm,  (B) 
2nm,  (C)  5nm,  and  (D)  llnm.  The  results  show  that  the  relaxation  is  much  longer  on  lnm 
particle  (~35ps).  The  relaxations  on  larger  particles  are  much  shorter  (~2.2ps).  The  fast  decays 
indicate  a  very  different  energy  relaxation  mechanism  on  the  larger  particles  - 
electron/vibrational  coupling. 
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Figure  6.  CO  vibrational  relaxation  dynamics  on  (A)  Pt,  (B)  Pd,  and  (C)  Au  nanoparticle 
surfaces  (step  surface  site).  (D)  CO/metal  interaction  strength  from  literature  (PRL.  1996,  2141). 


Figure  1.  (A)  Illustration  of  CO  surface  sites.  (B)  2D  IR  spectrum  of  the  CO  molecules  on 
various  sites  at  Ofs.  (C)  The  vibrational  relaxations  of  CO  on  various  sites. 


Figure  8.  Waiting  time  dependent  2D  IR  spectra  of  CO  molecules  on  (A)  Inm  Pt,  and  (B)  5nm  Pt 
surfaces  by  laser-exciting  CO  molecules  on  the  bridge  sites  and  detecting  the  responses  on  both 
step  and  bridge  sites.  The  cross  peaks  on  the  step  site  grow  much  slower  on  lnm  particle  than 
5nm  particle,  indicating  a  much  slower  heat  generation. 


7 


Surface  CO  Stretch 


>2nm 


Inm 


Electron- 

vibration 

coupling 


~50ps 


Vibratio  n  - 
vibration 
coupling 


Surface  I  Phonon 

Electron  (lattice  vibration) 

-  1  <2ps 

^  >10ps 
Environment 


Figure  9.  Energy  dissipation  pathways  on  Pt  nanoparticle  surfaces. 
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